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ABSTRACT 

We present the first measurements of the angular correlation function of galaxies selected 
in the far (1530 A) and near (2310 A) Ultraviolet from the GALEX survey fields overlapping 
SDSS DR5 in low galactic extinction regions. The area used covers 120 sqdeg {GALEX - MIS) 
down to magnitude AB = 22, yielding a total of 100,000 galaxies. The mean correlation length 
is ^ 3.7 ± 0.6 Mpc and no significant trend is seen for this value as a function of the limiting 
apparent magnitude or between the GALEX bands. This estimate is close to that found from 
samples of blue galaxies in the local universe selected in the visible, and similar to that derived 
at z ~ 3 for LBGs with similar rest frame selection criteria. This result supports models that 
predict anti-biasing of star forming galaxies at low redshift, and brings an additional clue to the 
downsizing of star formation at ^ < 1. 

Subject headings: Galaxies: UV - Correlation Function Evolution - Star Formation 
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1. Introduction 

In the current paradigm of structure forma- 
tion, the bulk of the most massive systems form 
in a cold dark matter-dominated universe by the 
merging of less massive units formed earlier. In 
parallel to this hierarchical evolution, recent ob- 
servations point to the so-called "downsizing", 
namely the fact that in galaxies having high bary- 
onic masses the bulk of stars formed at high red- 
shift (z > 1), while in galaxies having low bary- 
onic mas ses the bulk of s t ars formed at lower 
redshift (ICowie et all Il99d iHeayens et al.l [200I 
Bundv et al. 20061 : Jimenez et al. 2005 , and also 
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De Lucia etldl (|2006l ) and iNeistein et all (|2006l ) 

for results from simulations). The star formation 
efficiency shows a strong decline at < z < 1, as 
measured by the evolution of the star formation 
rate d ensity (^Hopkins fc BeacomI 



1996 



2000; 



ISch iminovich et al.l 
IWilson ( 



et alJ 



2006': 'Lilly el^ 



2005: Sullivan et al 



These epochs also 



see the bulk of the build-up of the bimodality 
in galaxy properties of the local univ erse, which is 
appar ent in their color distribution Baldrv et al] 



2004 . 



()2004l) . morphologies (fKauffmann et al 
spectral class ( Mada wick et al. 20021) and spatial 
distribution (jBudavari et al.l l2003h . Understand- 
ing the full picture is complex as this evolution 
is the res ult of the interpla y of several physical 
processes (|Faber et al.ll2005[) and combine the ef- 
fects of initial galaxy formation conditions ("na- 
ture") with galaxy evo lution events ("nurture") 
( Kauffmann et alJ I l2004h . [n this context, tracers 
that measure over cosmic time galaxy populations 
selected with homogeneous physical criteria are of 
primary interest. They help compare observations 
to simulation predictions over a large range of 
redshifts with reduced uncertainties, and allow a 
study of the redshift evolution of galaxy properties 
derived from different surveys. 

The ultraviolet (UV) range of the spectrum 
meets these conditions: UV luminosities provide 
a good mea sure of r ecent star formation within 
galaxies (Ke nniciifrO Il998l ). modulo attenuation 
by dust, and has been widely used at high red- 
shifts to study the properties of the Lyman Break 
Galaxies (LBGs) (Giavalisco fc: DickinsonI 2001 



Shaplev et cIll2003HSteidel et al.lll995^ . As large 



amounts of data are now becoming available at 
lower redshifts as part of the GALEX surveys 
( Martin et al. l2005h . the restframe UV spectral 
domain is presently well sampled over the full 
< z < 6 redshift range. Furthermore, compar- 
ison of results from high and low z UV-selected 
samples is eased by the fact that the UV luminos- 
ity dens ity fraction prob ed at high and low z are 
similar (jHeinis et al.ll2007l . hereafter Paper II), due 
to the strong lum inosity ev olution of the UV lumi- 
nosity function ( Arnouts e t al. 2005). Noticeably, 
during the epochs probed by GALEX the prop- 



^We define the UV luminosity density fraction of a given 
sample as the ratio of the UV luminosity density encom- 
passed by the sample over the total UV luminosity density 
at the same redshift. 



erties of active star forming galaxies show a very 
fast evolution. 

The wealth of UV-selected data now available 
at low redshifts enables statistical studies in the 
context of the downsizing of star formation, and 
in particular searches for links between the star 
formation properties and galaxy environment in 
terms of galaxy or dark matter density. Here we 
focus on the evolution with redshift of the link of 
star formation with Dark Matter and particularly 
the evolution of the class of Dark Matter halos 
hosting actively star forming galaxies since z ~ 1. 
This can be achieved by the study of the clus- 
tering of galaxies: at high redshift, LBGs stud- 
ies show that UV selec ted galaxies inhabi t high 



galaxy density regions ( Steidel et al.l Il998h . and 



are strongly biased with respect to the underly- 
ing Dark Matter, with more a ctively star form- 
ing galaxies being more biased (lAdelberger et al 



l2005MGiavalisco fc Dickinson" 200ll : lFoucaud et al 
i2003i . see Giavali sco (2002) for a review on the 
properties of LBGs). We propose to extend such 
studies to low redshifts using similar selection cri- 
teria. This paper is the first in a series and 
presents the methods and first results of angu- 
lar clustering measurements from GALEX data. 
The following section presents the datasets and 
the derivation of the redshift distributions. Sec- 
tion [3] presents two methods to derive the angular 
correlation function from a set of fields, and a dis- 
cussion about the behavior of these methods with 
respect to photometry inhomogeneity. In section 
|4] we present our results on the angular correla- 
tion functions and correlation lengths. To provide 
the crucia l link to Dark Matter halos, wc use the 
analytical IMo fc Whit3 (|2002) formalism that we 
present in sec. 14.31 We end with a short discussion 
in sec. [5] 

Throughout the paper a KCDM cosmology is 
assumed with matter density = 0.3, vacuum 
energy density 51 a = 0.7, and a Hubble parameter 
h = 0.7 where iJo = 70km s~^ Mpc~^ All corre- 
lation length values taken from the literature have 
been co nverted acco rdingly using equation (4) in 
iMagliocchetti et all (|200Q ). 

2. Primary fields selection 

We use Medium Imaging Survey (MIS) fields 
from the GALEX Release 2 (GR2), which al- 
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lows us to probe the clustering of faint sources 
in FUV and NUV at three limiting apparent AB 
magnitudes 22.0, 21.5 and 21.0. The magni- 
tudes we refer to are corrected for Galactic ex- 



tinction using Schlegel maps (ISchlegel et al.lll998l ) 
and the Cardelli et al. I (Il989l ) extinction law, un- 
less specified otherwise. The average color excess 
in the fields, derived from the lSchlegel et all maps, 
ranges from 0.08 to 0.12. The extinction coeffi- 
cients Afuv/E{B -V) and Anuv/E{B -V) are 
respectively 8.29 and 8.61. 

We start with 348 MIS GALEX fields over- 
lapping SDSS DR5 of which only a subset will 
be kept because of galactic extinction (see sec. 
13. 2p . We only include sources within 0.5 deg ra- 
dius from the field center, since artifacts concen- 
trate near the periphery of the field of view (see 
iMorrissev et al . 2005) and photometric accuracy 
decreases beyond this limit. We used only objects 
within the GALEX primary resolution. We used 
SDSS masks to mask out holes, bright stars and 
trails; we used also GALEX masks, as well as ad- 
ditional ones to mask out resolved galaxies or arti- 
facts not predicted by the GALEX pipeline. Only 
objects with an SDSS match within 4" are kept, 
and the closest match is used. Galaxies are as- 
sumed to be SDSS galaxies (type = 3 following 
the m orph ological classification o f iLupton et al 



J2OOI') and'Stoug hton et all (|2002l )). To check the 
effects of possible residual QSOs in our sample, 
we removed from our sample AGN-dominated ob- 
jects as objects classified QSOs by a template fit- 
ting procedure!!. Impacts on the results are found 
negligible, hence we do not remove QSOs classified 
objects from the sample in the following. 

2.1. Redshifts distributions 

To get the redshift distributions of the sam- 
ples for the different magnitude c uts, we use the 



polyn o mial fi t method described in lConnoUv et al 



(|l995l 119971) . This method requires first to be 



trained with a spectroscopic sample. We train on 6 
bands [NUV from GALEX , as all objects do not 
have FUV photometry, and the 5 SDSS bands) 
with 17,843 objects from the sample described in 
sec. [2] having SDSS spectroscopic redshifts. We 
then apply the coefficients derived from the train- 



ing set to the whole sample. 

We performed a simple correction for the broad- 
ening due to photometric redshifts errors by as- 
suming that the photometric redshifts errors fol- 
low a normal distribution independent of the ob- 
ject magnitudes and redshifts, with the standard 
deviation a = 0.03 measured using all the avail- 
able spectroscopic redshifts from the SDSS. We 
check that the standard deviation docs not vary 
with apparent magnitude using our photometric 
redshift estimation on GALEX fields with SDSS 
overlap and the independent and deeper spec- 
troscopy from IPapovich et al. ( 20061 ). 



Following lEfstathiou et al.l (|l99lf ). the parent 
distribution of the true redshifts is described by 
the following parametric shape: 



dN 
dz 



^A, 



exp 



(1) 



■^Le Phare; available and documented at this URL: 
|http : //www ■ oamp . f rZ sirnouts/LE.PHARE ■ html| 



We fit this shape convolved by a Gaussian with 
a = 0.03 to the observed photometric redshift dis- 
tribution. FigUshows the gaussian-convolved best 
estimate N{z), along with the measured distribu- 
tion for the 174 least extinct fields (see sec. 13. 2p . 
and table [1] lists the parameters of the true distri- 
butions. 

3. Choice of ACF estimation method 
3.1. Methods description 

Given that the GALEX data are extracted 
from relatively large numbers of similar expo- 
sures, the question arises as to the most appro- 
priate method to retrieve the available informa- 
tion. A first, straightforward approach to mea- 
sure the ACF from a group of non-overlapping 
fields is to treat all of them as disjoint subfields 
of one large, dis c ontin uous field and to apply the 
Landv fc Szalav ( 19931 ) (hereafter LS93) estimator 
on it: 

DD - 2DR + RR 

= RR 

where DD (resp. DR, RR) is the number of data- 
data (resp. data-random, random-random) pairs 
from all fields (including cross pairs from different 
subfields), normalized by the suitable pair num- 
bers. In the case of this composite field method 
(hereafter CF), the number of random points is 
fixed for the global field, and not for each individ- 
ual GALEX field. This is the ideal method which 
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in principle allows one to extract all the available 
information. In particular, this method reduces 
the integral constraint bias and the noise, espe- 
cially at large angular separations. 

Although best in the ideal case, the CF method 
requires precise homogeneity of the data and may 
not be robust in practice. We therefore introduce 
another estimator, which we define as the follow- 
ing pair-weighted average (hereafter PW) of the 
ACF measured in each field individually: 



wpw{0) 



j:^RR^{0)w,{e) 



(3) 



where Wi is the ACF estimated from field^ijilone 
computed with the LS93 estimator, and RRi the 
number of random-random pairs in the random 
catalogue constructed for this field (see Appendix 
[5] for a derivation of this formula (eq [3]) from 
wcf)- The RRi term involves pair numbers and 
field geometry information. The PW method is 
by construction insensitive to field-to-field fluctu- 
ations - and thus best suited for the peculiar MIS 
geometry. A drawback of the PW method is the 
increase of the integral constraint (IC) bias be- 
cause of the smaller angular extent of the fielcjfl as 
well as an increase in the noise. The integral con- 
straint can be relatively well corrected for using its 
estimate given by LS93. To compute it, we assume 
that the real correlation function is a power-law 
Aw9~^ and we fit A^O^^ - I{A^,5) to the data, 
where I{Ayj,5) = 1/9? j^-^Ayje^^dVlidQ,2, inte- 
grated over a GAL EX field. This method is simi- 
lar to that used bv lRoche fc EalesI (|l999l ). except 
that 5 is left as a free parameter. In the following, 
"PW method" will refer to the IC-corrected tech- 
nique. We have checked the accuracy of the above 
correction of the IC bias using a lO Odet?^ synthetic 
catalog derived fro m GallCS (Ha tton et al. 20031 



Blaizot et all 120051 ). The ACF has been com- 



puted with the CF and the PW methods from 
50 randomly positioned fields of radius 0.5°. The 
GallCS-magnitude cut was chosen to obtain ap- 
proximately the same mean number of galaxies as 
found using the NUV < 22 cut. The results of the 
CF and PW methods have been found undistin- 
guishable for the model catalogs (fig [2]). 



•^In the case of independent fields, the IC in the PW method 
is typically higher than that of the CF method by a factor 
of the number of fields. 
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Fig. 1. — Derived redshift distributions for the 
two GALEX bands using three magnitude cuts (his- 
tograms). The solid curves show the best fitting N(z) 
convolved by a Gaussian with a — 0.03 (see text). 
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Fig. 2. — Validation of the method used to correct for 
the Integral Constraint bias. The ACF is computed 
from 50 randomly chosen fields in a synthetic catalog 
with the CF (filled squares) and the PW (filled circles) 
methods. The dashed (resp. solid) line shows the best- 
fit of the PW result uncorrected (resp. corrected) for 
the Integral Constraint bias (see text for details). 
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Fig. 3. — Angular correlation function for the CF (left) and PW (right) methods (see text). Upper panels show 
FUV < 21.5 selection and lower panels NUV < 21.5. The solid line shows the ACF for the whole sample; the other 
curves represent the ACF of the fields according to the probability that their magnitude distribution is drawn from 
the sample distribution than the whole sample (p < 0.01, dot-dashed; p < 0.05, dashed; p > 0.05, dotted.) 
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3.2. Systematic effects 

To test the sensitivity of the CF and PW meth- 
ods to systematics, we used a statistical approach 
to decide whether the photometry of a given field 
is drawn from the same distribution than the pho- 
tometry of the whole sample. To this aim, we used 
the Mann- Whitney test, which is independent of 
the size of the input samples; we use only objects 
with type= 3 during this process to avoid strong 
star contamination. For each field, we create a 
test sample built from the magnitudes of the ob- 
jects in this field, and a control sample from the 
magnitudes of the objects that belong to all other 
fields. The Mann- Whitney test provides as out- 
put the probability that these two distributions 
are the same. We show figure |3] the ACFs of fields 
grouped according to their value of this probabil- 
ity, using the CF method (left), or the PW method 
(right), for FUV < 21.5 (top) or NUV < 21.5 
(bottom). The results obtained from the CF and 
PW methods show significant differences. The CF 
method results show an excess of power especially 
at large scales; the amplitude of this excess in- 
creases as the probability that the field photome- 
try is the same than the overall sample photom- 
etry decreases. Conversely the results of the PW 
method are fairly insensitive to photometry inho- 
mogeneities; there is an overall power excess at 
FUV < 21.5 for the fields with a probability lower 
than 0.01, but the ACF of the whole sample is very 
similar to the ACF of the best fields (p > 0.05). 

There are several sources of systematic errors, 
which, although their individual effects are weak, 
may, combined with each other, yield the trends 
observed. Similar trends are observed when the 
fields are binned according to the mean Galac- 
tic extinction. However, the cross-correlation be- 
tween galaxies and dust maps using both CF 
and PW methods is found at least 5 times lower 
than the autocorrelation at scales where the lat- 
ter is positive {6 < 0.2 deg), and no obvious trend 
was found between the amplitudes of this cross- 
correlation and the Galactic extinction. On the 
other hand, the amplitude of the cross-correlation 
function between galaxies and background maps is 
higher in fields with higher mean Galactic extinc- 
tion. Inhomogeneities may also arise from pho- 
tometry drift with ti me, but sourc es drifted less 
than 0.1 magnitude llMorrissev et al. 2007) over 
the whole GALEX mission; a drift of this am- 



plitude has a small effect on the CF method, as 
expected from tests on mock catalogs. Studying 
sources observed several times in overlapping re- 
gions shows that field-to-field fluctuations are less 
than 10 % beyond what is expected from Poisson 
statistics. Note however that this result is based 
on a few sources per field located at the edges of 
the field, where photometry accuracy decreases. 
Star contamination can lower the amplitude of the 
PW method, as an addition of an uncorrelated 
population, while it may contribute to the effects 
observed with the CF method, given the variations 
of star counts with Galactic latitude. According 
to template-fitting based classification, the frac- 
tion of stars in SDSS objects with type = 3 is 2 
% in NUV and 8 % in FUV; we checked that this 
has small effect on the PW method. 

All these tests suggest that there is some source 
of field-to-field variations in our data, likely due 
to a combination of zero-point calibrations, back- 
ground fiuctuations (correlated with Galactic ex- 
tinction), etc ... The PW method is fairly insen- 
sitive to any systematics, as expected, and we are 
thus confident that it is a robust estimator. It 
is this method we chose to use in the rest of the 
paper. 

Conservatively, for the remainder of the paper, 
we restrict the analysis to the 174 fields with the 
lowest Galactic extinction {{E{B - V) < 0.04}). 
The number of galaxies at the different limiting 
magnitude cuts are given in table [TJ 

The characteristics of the UV dust attenua- 
tion in galaxies are not known to be correlated 
with the large scale structure or the galaxy den- 
sity (even in the extreme cases of clusters, see 
e.g. iBoselh fc Gavazzil |2006|) : the effect of inter- 
nal dust has thus been taken as an uncorrelated 
noise source on the UV fluxes and its effect on 
the ACF neglected. This allows direct comparison 
with clustering studies of high redshift restframe 
UV-selected galaxies. 

4. Angular correlation function and corre- 
lation length 

4.1. Measurements 

We use the PW method described in sec. [3] to 
measure uj{9) from the 174 fields with the low- 
est Galactic extinction, using logarithmic-width 
bins of A log 6* = 0.27, and 6'™„ = 0.005° and 
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Table 1 

Sample description, power-law best fits parameters and comoving correlation lengths 







FUV 






NUV 




Limiting mag 


22. 


21.5 


21. 


22. 


21.5 


21. 



An, X 103 [dcg*] 
S 

ro [Mpc] 



44651 
0.18 
5485.96 
0.16 
2.18 

0.80±0.05 
4.2^ 



22655 
0.15 
3616.28 
0.13 
2.29 

q 7+3.7 
a. (2.7 

0.80±0.08 



11418 
0.12 
2417.72 
0.11 
2.76 

0.75±0.12 



1+0.5 
-0.4 



3.7 



+0.7 
-0.6 



2.8 



+0.9 
0.7 



99368 

0.25 
8255.5 

0.18 

1.67 

q 9+0.5 
■^•^-0.5 

0.89±0.04 
4.0^ 



^+0.3 
-0.3 



48274 

0.21 
5084.1 

0.13 

1.61 

0.88±0.05 
3.7"^ 



7+0.4 

-0.3 



22948 
0.16 
3291.14 
0.1 
1.66 
6 3+° '' 
0.84±0.09 
3.3"* 



3+0.7 
-0.5 



* Number of galaxies in the samples 
tMean photometric redshift 

t Parameters of the true best fit redshift distribution (see text) 

Note. — The amplitude and slope of best fit power laws to the angular correlation function, and 
hence the comoving correlation length account for the Integral Constraint correction (see text). No 
attempt to remove residual QSOs from photometric rcdshifts is performed here. 



&max = 0.4°, which probes scales in the range 
0.05 to 4 comoving Mpc at the median redshift 
{z — 0.15) of the samples considered here. The 
results are plotted in fig.|H the la errorbars repre- 
sent internal scatter derived from jackknife resam- 
pling of the 174 GALEX fields used for the ACF. 
In order to check for any instrumental contribu- 
tion to the ACF such as residual non uniformities 
of the sensitivity across the field of view, the PW 
method has been applied to stars, selected as ob- 
jects whose SDSS counterparts with type = 6. For 
stars we find no significant deviation from a null 
correlation function. 

We fit the results using the method described 
in sec. [3l Our best fits for the different samples 
are given in table [U where 9 is expressed in de- 
grees. The error bars on and 6 are the pro- 
jected {xliin + 1) contour. 

To derive the comoving co rrelation lengt h, rg, 
we used the Limber equation (|Peebledll980h with 
the true deconvolved redshift distributions (see in 
sec. I2.1|) . The results are given in table [TJ The 
uncertainties on ro have been assumed to be the 
extreme excursions of tq in the projection in the 
(ro, S) and (rp, planes of the contour at the 



68% probability in the {A^, S) plane (see figE]). 

4.2. Comparison with previous studies 

Given the error bars, the slopes 5 found 
for the different magnitude cuts in the two 
bands are compatible with a constant value 
S ~ 0.81 ± 0.07. This is steeper than re- 
ported i n several stud ies base d on blue galaxies 
at low z (Budavari et ah, 200^: Zehavi et al. 20021 : 



Madgwick et al.ll2003f ). and restf rame UV-selected 
galaxi es at higher redshifts (|A - delberger et al 



20051 IPorciani fc Giavaliscol 1200^, all of them 
consistent with a value oi 5 ~ 0.6. How- 
ever our meas ur emen t is in agreement with 
Iciavahsco et al. ( 19981): Giavahsco fc Dickinsonl 



( 2001 ): lFoucaud et all (|2003l ). Moreover JCoil etaP 
(j2004[ ) noticed a steepening of the slope not only 
for the reddest but also for the bluest galaxies of 
their samples. We discuss in Paper II the depen- 
dence of 6 on UV luminosity. 

With an average comoving correlation length 
3.7±0.6 Mpc at z - 0.2, the present GALEX data 
sets confirm the low clustering of the rest-UV se- 
lected galaxies in the local universe observed by 
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Fig. 4. — Angular correlation function measured in the 174 GALEX fields with the lowest Galactic extinction, for 
the FUV (left) and NUV (right) GALEX bands at three magnitude cuts: circles, muv < 22; squares, muv < 21.5 
and triangles, muv < 21. Dashed lines show the power law best-fit uncorrected for Integral Constraint. The upper 
axis shows the comoving distances corresponding to the angular scales ai z — 0.15. No attempt to remove residual 
active nuclei by photometric redshift template fitting has been made beyond the SDSS classification. 



Heinis et al.l (|2004l ). The new mean value is 25% 
lower, though both measurements agree within er- 
ror bars. Assuming the average values of ro and 
S quoted above, the corresponding bia s defined at 
8 Mp c by bs = o^.gja^ 



(e.g. iMagliocchetti et al 



0.2, a significant anti- 



2000) is 0.61 ± 0.09 at z 
bias. 

It is well-known that blue galaxies are less 
strongly correlated than red ones, and not sur- 
prisingly the small correlation length found in this 
study is comparable to that measured for blue 



( 20041 ) 



galaxies in the local universe: iCoil et al. 
report a comoving tq of 2.54 ± 0.37 Mpc for the 
class of blue galaxies defined by (0.2 < R — I < 
0.4) in their visible-selected sample, which spans 
the redshift range 0.3 — 0.6, with the lowest cor- 
relation length among all of their galaxy sub- 
samples. The GALEX restframe UV selected 
galaxies are nevertheless even less correlated than 



the galaxy cla s s T4 (bluest of 4 classes) from 
Budavari et al.l (120031) for which the y deriv e a ro 
of 6.44 ± 0.27 Mpc. iHawkins et"al] (|200lh com- 
puted the redshift-space correlation function from 
far-infrared selected galaxies in the local Universe 
{z ~ 0.03). Converted to real-space, their esti- 
mate of the correlation length (rg = 5 ±0.33 Mpc) 
of the hotter galaxies, i.e. the most star-forming, 
is higher than ours from UV-selected galaxies. 

Low-redshift restframe UV-selected galaxies 
possess correlation lengths slightly lower than 
those derived from high-z restframe UV-selected 
samples (see fig. [6|). Note that the comparison 
with results from higher z samples is not straight- 



(IGiavalisco & Dickinson! 


2001 : lAdelberaer et al. 


2005: Zehavi et al. 2005 


: Norbere et all 2002h: 



brighter objects are more clustered than fainter 
ones. The GALEX samples are the faintest of the 
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Halo masses [Log [M^JM^] 



10.0 




Adelberger et al (2004) □ 
Arnouts et al (2002) * 
Foucaud et al (2003) 
Giavallsco & Dickinson (2001) o 
Helnis et al (2004) o 
GALEX FUV • 
GALEX NUV ■ 



0.001 



0.010 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



Fig. 5.— x~ 

on S and m for the NUV < 22 result. Top left: 
contours of constant in the {Aw,d) plane. The inner 
(resp. outer) solid line corresponds to the 68.3% (resp. 
95.4%) confidence level. The dashed line shows the 
Xrnin + 1 coutour; its projcctious on the axes give the 
error bars on Aw and 5. Top right: comoving ro as 
a function of 5 given using the values of y4„ and 
6 included in the 68.3% confidence level. Bottom left: 
comoving ro as a function of A^ given 5 using the 
same values. The errors bars on ro are the extrema of 
this distribution. The solid line intersections show the 
location of the best fit in each panel. 



restframe UV-selected samples considered here: 
the mean absolute magnitudes of the FUV and 
NUV samples are Mpuv — —18.3 and M^uv = 
— 18.8, while the Lyman b reak galaxies sample s 
of lAdelberger et al.l (l2005l): lArnouts et all (|2002|): 



Foucaud et al.l (|2003[ ): iGiavalisco &: Dickinson 
(l200l|) are all brighter than Muv = -20. We 
study in details the luminosity dependence of clus- 
tering within the GALEX samples in Paper II. 

4.3. Comparison with Dark Matter Halo 
clustering predictions 

In this section we use the formalism described 



by iMo k, White! (|2002l ) to compute the correla- 
tion length of Dark Matter Halos (DMHs) above 
a given mass as a function of redshift. We assume 
that the spatial correlation function of DMHs with 



Fig. 6. — Comparison of the evolution with redshift 
of the correlation lengths obtained from restframe UV 
selected samples with the correlation lengths of Dark 
Matter Halos more massive than Mmin (color-coded). 
The inset shows the low z points. 



masses greater than Mmin at a redshift z is well 
approximated by a power-law: 



-Ih 



(4) 



where ro,SMm.in., z) is the cor r elatio n length of 
such halos. The IMo fc Whitl (|2002l ) formalism 
provides analytical equations for the abundance 
and the bias factor of the halos, n(M, z) and 
b{M,z), given their mass and redshift. The ef- 
fective bias of the halos more massive than a min- 
imum mass Mmin at a redshift z is then given by: 



heff{Mmin,z) 



J^^^^ h{M, z)n{M, z)dM 
/aT , «(M,z)dM 



(5) 



The rms density fluctuations of the halos is linked 
to the rms density fluctuations of the underlying 
mass at 8 /i~^Mpc by: 

Crs.,h{Mmin, z) = beff{Mmin, z)a{z)s.m 

where the subscripts h and m denote respec- 
tively halos and underlying mass, a nd <j{z)s.m = 
a{0)s,mD{z) (see lMo fc Whitell2002l ) withcrfO)« ^ = 
0.9. The correlation length of the DMHs with 
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mas ses M > Muin at z is t hen obtained using 
(e.g. iMagliocchetti et aLllioOol ): 



expecte d to be stronger for st ar forming galaxies 



<yS,,h{Mmin, z) = 1 / C 



7h 



(Mmin, Z) 



(6) 



where = 72/[(3 - 7)(4 - 7)(6 - 7)2^]. We 
assumed that the slope of the spatial correlation 
function of the halos is 7^1 = 1.8, after having 
checked that the results are rather insensitive to 
the adopted value if 1.5 < "fh < 2.5. Fig. [6] shows 
the redshift evolution of the correlation length of 
DMH with masses lO'^ < M™„ < W^Mq. 

In the framework of Hal o Occupation Distri 



bution (HOD) models (e.g. iBerHnd fc Weinberg 



2002t ICoorav fc ShethI I2OO2I ) . recent studies have 



pointed out that the galaxy correlation function 
is likely to be the sum of two components. The 
first component dominates at small scales, de- 
scribing the correlation of galaxies that are in 
the same halo, and the second component ac- 
counts for galaxies in different halos, dominating 
at large scales. The present sample however does 
not provide sufficient constraint to fit HOD mod- 
els to our results — we plan to perform this in fu- 
ture papers. Assuming that the correlation func- 
tion of DMHs is a power-law at all scales and 
that each halo hosts at most one galaxy, a di- 
rect comparison of the GALEX results with the 
correlation lengths of dark matter halos (fig O 
shows that the UV selected galaxies in our sam- 
ples have the same correlation lengths as ha- 
los with masses lower than Mmin — lO^^M©. 
At z > 2, clustering measurements from LBGs 
samples show that halos with comparable clus- 
tering strengths ha ve Mmin ^, lO^^Mc^, a s al- 
ready m entioned by Adelberge r et al. ( 20051 ) and 
[Ciavalisco fc Dickinson (2001). These results sug- 
gest that the characteristic mass of halos hosting 
active star formation has decreased from z — 3. 
Note that taking luminosity evolution into account 
does not weaken this result, since at low redshifts 
UV-selected samples actually probe the same UV 
LD fraction as their high redshift counterparts (see 
Paper H) . This mass evolution, as well as the bias 
evolution can be interpreted as additional evidence 
for the "downsizing" scenario of star formation, 
though it applies here to the mass of the underly- 
ing halo rather than to the baryonic mass. Note 
that t hese is a correlation b etween halo and galaxy 
mass ( Shankar et al. 200d ). though its scatter is 



(see e.g. lYoshikawa et al.ll2001l ). Theoretical stud- 
ies also predict tha t the SFR increases with halo 
mass (e.g. at z = S. lBouche et al.J 20051): however 



in prese nce of AGN feedb ack(e.g. lDi Matteo et al 
[2005; Croton et al.ll2006l ). or when taking into ac- 



count gravitational heating (jKhochfar fc Ostriker 
20071 ). this trend reverses at lower redshift. 



The conclusions of this work are developed in 
Paper II. 

5. Conclusions 

We presented here the first clustering measure- 
ments from the GALEX data. These data provide 
a unique basis to statistical studies of star for- 
mation in galaxies at low redshift from their UV 
continuum. The same tracer can now be used in 
an homogeneous way over a large redshift range 
(0 < z < 4) to investigate the processes driving 
star formation evolution. We discussed the impact 
of the Galactic extinction on the clustering mea- 
surements and used a method insensitive to the 
color excess. We measured the clustering by the 
angular correlation function, and fitted our results 
with a power law parametrisation: w{9) = AiijO~^ . 
We derive steep slopes, S ~ 0.81 ±0.07. Assuming 
photometric redshift estimation, we compute the 
correlation length, tq . The results confirm the low 
clustering of UV-selected galaxies at low redshift 
(ro = 3.7 ±0.6 Mpc). Comparison with analytical 
modeling shows that active star forming at z < 0.4 
present the same correlation lengths than DMHs 
with Mmin < 10"'^'^ M0. This result is in agreement 
with the "downsizing" scenario. 

It is with great pleasure that we thank Jean- 
Michel Deharveng for support and discussions. 
GALEX (Galaxy Evolution Explorer) is a NASA 
Small Explorer, launched in April 2003. We 
gratefully acknowledge NASA's support for con- 
struction, operation, and science analysis for the 
GALEX mission, developed in cooperation with 
the Centre National d'Etudes Spatiales of France 
and the Korean Ministry of Science and Technol- 
ogy- 
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Pair- Weighted Average estimator for the Angular Correlation Function 

In this section we discuss the weighted estimator presented in sec. [3l The expression can be derived 
directly from the definition of the LS93 estimator with some assumptions. The LS93 estimator is: 

DD{e)-2DR{e)+RR{e) 

= mm 

where DD, DR and RR are normahzed by the suitable pairs number : 

UgiUg - 1) 

DR ^ ^ (3) 

RR = (4) 

where Ug is the number of galaxies in the sample and is the number of random objects in the random 
sample. In the following we do not recall the 6 dependence of the different quantities. 

Let us consider the case of the CF method (see sec. [3]) applied on N fields positioned on the sky in such 
a way that no cross pair between objects from different fields has to be accounted for in the computation 
of w{9). The total number of pairs over all the fields in each angular bin can then be expressed using the 
number of pairs in each field: 

i=l i=l 

where DDi is the number of data-data pairs and Ug^ the number of galaxies in the ith field. The same 
equations hold for DR and RR. 

When computing the ACF of one field individually, we consider 100 random samples with the same 
number of random points that galaxies in this fielcQ. RRi is then the average of the 100 computations. So 



Ug. = — rii and Ug — Ur — n. Then the LS93 estimator can be written: 



w 



2RB. 



^ {DD + RR) - \dR 



n{n — 1) 



n(n— 1) 

Let us consider the term in brackets; with our assumptions it yields: 



(6) 



n{n — 1) 

then introduce the term RRi/ RRi in both sums: 



2 v- ^'K^ ^\ dd, + rr,)-^Y.<^^ (^) 



n{n — 1) 



In the case of the CF method, the total number of random points would also be fixed to rig, but the number of random points 
in each field is allowed to be different of Uo ■ 
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The ACF of the ith field can be written as: 



Hence [5] becomes 

2 n^iui - 1) 



wi- + + W3. with 



DRi 



RRi 
RRi 



RRi 



n{n — 1) 

At this stage we also assume that 1 so that ni{ni — 1) ~ ri|, and hence n(n — 1) ~ n^; [5] yields 



Coming back to eq. [5] we finally get 



^YrR^w, (10) 



Wpw(O) = — 11) 
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